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ABSTRACT 



^ The nonlinear oscillator model is useful to basically understand the most important properties of 

in 

O nonlinear optical processes. It has been shown to give the correct asymptotic behaviour and to 

O provide the general features of harmonic generation to all orders, in particular dispersion relations 

o^ 

^ and sum rules. We investigate the properties of Pump and Probe processes using this model, and 
^ O study those cases where general theorems based on the holomorphic character of the Kubo response 
>^ functions cannot be applied. We show that it is possible to derive new sum rules and new Kramers- 
Ph Kronig relations for the two lowest momenta of the real and of the imaginary part of the third order 
susceptibility and that new specific contributions become relevant as the intensity of the probe 
increases. Since the analitic properties of the susceptibility functions depend only upon the time 
causality of the system we are confident that these results are not model dependent and therefore 
have a general validity, provided one substitutes to the equilibrium values of the potential 
derivatives the density matrix expectation values of the corresponding operators. 
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I. INTRODUCTION 

The theoretical and experimental investigation of systems responding to Pump and Probe 
optical beams is an important branch of research in solid state physics [1,2,3,4,5,6]. The most 
important nonlinear effect is the modification of the optical response near each resonance, which is 
often called the dynamical Stark effect [7,8,9,10,11]. Other new effects are however present, such as 
the two-photons absorption [12,13,14,15] and the stimulated Raman scattering [16,17,18,19,20]. To 
completely interpret the above phenomena detailed calculations are necessary [21,22,23,24,25] but 
in general they are very difficult to perform with the required accuracy for realistic systems. For this 
reason general properties of the nonlinear response functions are very useful for the interpretation of 
experimental data and the testing of appropriate models [26,27,28]. hi particular, sum rules and 
Kramers-Kronig (K.K.) relations are expected to be of great help also in nonlinear optics, because 
they only depend on general properties, such as time causality, and must be verified by any physical 
system. 

Sum rules and K.K. relations have been derived for those specific nonlinear processes 
represented by the susceptibility z " {<^i,<^2, -(02,--'^2, ~<^2)' which is the most relevant for weak 
probes (whose frequency is (o) since it is proportional to the {n - 1)/2* power of the intensity of the 
pump beam. In fact it has the property of being holomorphic in the upper complex half plane of the 
variable cOj thanks to Scandolo's theorem [29]. 

When we consider nonlinear effects also due to the probe beam the previous treatment is not 
valid anymore and new poles in the n* order susceptibility functions appear in the upper complex 
plane. Scandolo's theorem in this case is not verified so that nothing can be said about K.K. 
relations and sum rules from the previous analyses. This is clearly explained in the book by 
Peiponen et al. [30] 

The main purpose of the present work is to study in detail the third order susceptibility 

Z^^^i^l) with a probe beam and a pump beam, and to derive general properties such as K.K. 

relations and sum rules. Since such general properties do not depend on the microscopic decription 
of the nonlinear system, we adopt a simple anharmonic oscillator model along the line described by 
Bloembergen [1] and Garrett [31,32] and previously adopted for the studty of harmonic generation 
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[33] and for general studies of sum rules in special cases [34,35,36]. Such a model allows a 
derivation of the most important features, considering all nonlinearities and including also the 
nonlinear effects arising from the probe beam. We will show that new K.K. relations and modified 
sum rules can be obtained, taking into account optical rectification and considering the poles in the 
upper complex plane, which originate from the probe nonlinearities. 

From previous experience with harmonic generation processes [33] we feel confident that 
the results here obtained can be extended to the quantum mechanical analysis by substituting the 
expectation values of the potential derivatives to their values at equilibrium used in the classical 
description. 

In section n we describe our model and give formulas for x^^^ (^i )> while the higher orders 
expressions are given in the Appendix. We analyze the structure of the third order susceptibility and 
we show that the holomorphic and non holomorphic contributions can be separated and their 
respective properties can be considered. In section HI. we derive our modified sum rules. In section 
IV we derive the K.K. relations and discuss their connection with the sum rules. In section V we 
present our conclusions. 



II. THE NONLINEAR OSCILLATOR MODEL AND THE THIRD ORDER 
SUSCEPTIBILITY WITH NON HOLOIMORPHIC CONTRIBUTIONS 

We adopt the Lorentz oscillator model with a general potential energy mV(x) containing 
anharmonic terms, introduce a damping y and consider two incident optical beams, one used as a 
probe of frequency cOj, and one as a pump beam of frequency : 

E{t) = E^e-''"^' +E2e-''"^' +C.C. , (1) 



where c.c. indicates the complex conjugates of the previous terms. 

Expanding the potential around the equilibrium position we obtain the following equation 
of motion: 
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x + yx + colx + ^ 



«=3 



ok" 



X"-' eEit) 



{n-\)\ m 
^ ^ 



The solution of this equation can be obtained with an iterative procedure as described in ref. [1,2] 
and recently reported in detail in the book by Peiponen et al. [30]. 
We obtain for the linear contribution: 

P'''{co^)^E.x'''{coy j=l,2 (3a) 
with the usual : 



Z^\o)j) = —^ (3b) 



where : 

d((Oj )= 0)1-0)]- iro)j . (3c) 

A general expression of the nonlinear Polarization F^"^{rQ)^ + 30)2) is given in the 
Appendix A, where also a formal analysis of the properties of such functions in the complex cOj 
plane is given. We here concentrate on the third order susceptibility, and wish to prove in this 

particular case that general results like K.K. relations and sum rules can be extended to include the 
contributions which are nonholomorphic in the upper complex cOj plane, provided the poles are 
explicitely considered. 

The expression of the Z^^^^^i) can be obtained in the usual way by successive iterations 
[1,30,33] and also follows from the general result given in the Appendix A substituting n=3, r=l 
and s=0 in (Al). We obtain: 

j*^' («i )=z^^^ («i ; «i -«2 ' ^2 )+ z^^^ («i ; «i >o(- «2 ' ^2 ))+ 

+ Z^^^{co,;Q)2,o),-co2)+z^^^ {(^i ' + ^^2 ) 
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+ x^'^^ {(^1 \(Oy,(o^ -(o^ )+ Z^^^ (^^1 ; Oil '0(fyi -(o^ ))+ x^'^^ (oJi -(o^ ,2(o^ ) 



(4a) 



where 



X^^\o),;co,-co^,co^)=- 

X^^\co,;±co^,(o, +(o^) = 
X^'^{(0,;q),-q)„(0,)=-^ 
X^'KiOi--iOifi{(o,-o)J)= 
X^^^{co,-(jo„2(jo,)=^ 



d'Vix) 
dx' 

d^V{x) 



\ESe'Nlm' 



^ D(Q)y D(co,)D(-co,) 



\Efe'N/m' 



^Q)lD(a),fD(Q),)D(-Q),) 



\E}^e'Nlm' 



\Efe'N/m^ 



dx" \^(olD{(0^)' D{-(0^) 



d^V{x) 



\E,\^ e^Nlm^ 



2[ dx^ \^ D{(o,y D{-(o,)D{2(0,) 



(4b) 



(4c) 



(4d-4e) 



(4f) 



(4g) 



(4h) 



We observe that both holomorphic contributions, which obey Scandolo's theorem [26,29], 
and nonholomorphic contributions are present. Among the holomorphic contributions, the term (4b) 
gives the usual correction to the linear resonance (dynamical Stark effect), the term (4c) originates 
from optical rectification and is usually neglected, while the contributions (4d) and (4e) include the 
two-photon absorption and the stimulated Raman scattering. We notice that the holomorphic terms 
(4b) and (4c) together constitute the main contribution x^^^ ^\ )> holomorphic in the upper complex 
plane, as described in the Appendix A. The above described contributions have been already given 
in the literature [2,30] and are discussed in ref. [30], where also the problems due to the 
meromorphic terms (4f), (4g) and (4h) are enphasized. We observe that the nonholomorphic terms 
are due to the nonlinear effects of the only probe beam, (4f) corresponding to frequencies mixing, 
(4g) to optical rectification and (4h) to the two-photon absorption of the probe beam. 

Since the last three terms are nonholomorphic in the upper complex cOj plane their poles have 
to be explicitely considered in order to derive sum rules and K.K. relations for the total X^^ {o)^\ 
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We notice that two types of contributions appear in the holomorphic terms and in those which are 

not holomorphic, one proportional to |<?^F(x)/(?'^xJ^ , the other to [3'^F(x)/(?'^x]^ . The first is due 

to the sum of three frequencies, the second to the sum of two frequencies, one of which is already a 
combination of two. These are the only way possible ways to obtain a third order susceptibility 
Obviously, if the system has particular symmetries, some terms can be equal to zero; for 
instance in the case of inversion symmetry all the odd derivatives of the potential are zero so that 
only terms (4b) and (4f) survive. 

We can notice that the only poles of the upper complex plane appearing in the last three 
terms (4f), (4g) and (4h) correspond to the solutions of D{- ftJ^ ) = , i.e. : 

co,=±c + id = ±{q)1 - f + ir/2 . (5) 

We consider the above poles in order to obtain sum rules and K.K. relations of general validity for 
For convenience we present in fig.l and fig.2 a visual description of the terms 

2 2 

proportional to IE2I and lEil respectively. 



III. SUM RULES 

The general purpose of sum rules is to find the values of the moments of the real and of the 
imaginary part of the susceptibility. This is possible up to a given moment, above which the integral 
diverges. Before deriving the sum rules by direct calculations we want to show how it is possible to 
find linear relations between the four sum rules referred to the first two moments of the third order 
susceptibility z^^^ (^1 ) just by considering the following fundamental property of all susceptibility 
functions: 

Z'"i-o)*)={z'"i(o)j (6) 

This property requires that if Z^"\^) has a pole in a, it must have another pole in -a*, and that 
IS even while Imix'^"^ (q))) is odd with respect to co if co belongs to the real space. So we 
obtain for the principal parts of the integrals on the real © axis: 
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+00 +00 

P ^dQ)(0^'"Z^"^{a))=2 dQ)Q)^"'Re{x^"^{Q))) , (7a) 

— oo 

oo +00 

P ^d(0(0^'"^'z^"^ {(o)= 2i P j d«ft)""^'lm(;ir<"' {(o) ) , (7b) 

— oo 

for any integer m allowed by the convergence of the integrals. 

Now we restrict our analisys to the case n=3, where some conclusions can be easily 
obtained. 

Considering our expression (3c), we observe that D{-co*) = D{co )* since it is proportional to the 
inverse of a first order susceptibility (3b). As a consequence we find that the two residues of the 
poles (5) are one the opposite of the complex of the other. We denote the values of the two residues 
as a+ib and -a+ib , corresponding respectively to the poles at c+id and -c+id. We define So and S2 
respectively the zeroth and the second moment of the real part, while Si and S3 the first and the third 
moment of the imaginary part of the x'iH (^1 )' i-^- nonholomorphic contribution to the total 
susceptibility X^^^{^i) • Taking into account that the only the non-zero value of the moments of the 
holomorphic contributions is the third moment of the imaginary part [26], when we consider So, Si 
and S2 we can write x^^^ (^1 ) instead of x'^nl. (^1 ) without any difference. We then obtain: 

p"\dQ),x^'\Q),)=2 P\RQ(x^^\Q))=2S^=2m{a + ib-a + ih)=-A7tb , (8) 

—00 

P ]dco,coa''' {co,)= 2i P j dco,co,lm(x''' («i ) )= 2/*5i = 
— 00 ■ ^ * 

= 2m{a + ib\c + id )+ 2m{^ a + ib)^c + id ) = A7ti{ac - bd) 

We have found how to express the real and imaginary part of the values of the residues of the 
susceptibility function as linear combinations of So and Si. Since S2 and S3 are functions of the 
residues values and of the poles only, it is clear that it is possible to obtain for them an expression 
which depends only on the two lowest moments. The explicit calculation for S2 gives: 
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P \dco^O)fz^^^ {coJ=2 P\ dft)ift)f Re(;ir^^^ {co, ))= 2S^ = 2m{a + ib\c + idf + 
-oo • (10) 

+ 2m{- a + ib%- c + idf = -A7t{2acd + b{c^ - d^ )) 

Substituting the values of a and b as obtained in (8) and (9) we obtain: 

S^=So(c^ +d^)-2dS, (11) 

which, using the expression (5) for c and d, becomes: 

S,=a)X-}^, (12) 

Performing a similar calculation for the third moment of the imaginary part we obtain for the non- 
holomorphic contribution zi^li^i)'- 



P \dco,colx^l {co, )= 2i P I (oflm{x!^l {(o, ))= 2/^3= 2m{a + ib\c + idj + 

—00 

+ 2m{- a + ibX- c + idf = 4m{ac{c^ -3d^)- bdi^c^ -d^)) 



(13) 



Performing the same substitutions as above we obtain: 

S^=2Sod(c^ +d^)+S,(c^ -3d^) (14) 
Using expression (5) we have: 

S,=S,rQ)',+sXco'o-f) (15) 

Having established the above general results we proceed to the direct calculations of the 
moments of the third order susceptibility. We observe that we can express the sum rules as 

p]do)iO)f"'^' Imix'^'^ {o)^ )) (16) 
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and 



(17) 







with m=0,l since higher moments of the susceptibility would diverge because the terms (4b) and 
(4c) asymptotically decrease as co^"^ . We perform the integrals of the terms (4f), (4g) and (4h) 
closing the contour in the upper complex cOj half-plane and summing all the residues of the poles 
inside the semicircle considered, while we use the K.K. relations and the superconvergence theorem 
to compute the contributions of the holomorphic terms [26,27], obtaining immediately that all the 
moments of the terms (4d) and (4e) are zero because of their fast asymptotic decrease (Q)^''). After 
lenghty calculations we finally obtain the general sum rules of the total susceptibility ;jf {co^ ) 








(18) 




\E,f {l4 f +5coiy N/m' \ d'V(x) 



(19) 



4Sy'Q)o(2f +0)1) [ ok' 




+ 



n\ 



(20) 
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7r\E^^ e'^ Njnr' 



(21) 



We underline that the only parts proportional to the intensity of the pump beam are the first two 
terms of the third momentum of the imaginary part, which originate from the holomorphic 
contribution J*^^(<Wi). They correspond to those previously derived [27]. However another term, 
proportional to the probe intensity, appears in this moment. This third term is S3 , while we can 
identify (18) as So , (19) as S2 , (20) as Si. So all the Si with i =0,1,2,3 are proportional to the probe 
intensity. We have also the confirmation that all these terms obey the linear relations (12) and (15), 
as can easily be checked with algebraic calculations. 



IV. KRAMERS-KRONIG RELATIONS 



It is useful to verify the existence of specific K.K. relations and to show how they are 
modified with respect to those already obtained for the holomorphic contribution x'^^^ (^1) [26]. 

The method is to consider explicitely the poles in the upper complex plane and to perform 
the Cauchy integrals for all the moments which are allowed by the asympotic behaviour. The 
asymptotic behaviour of is ^ Q)^ , so that the moments to be considered are ^iid 

Q)fx^^\cji)^). We write explicitely the two Cauchy integrals, considering the poles at (c+id) and at (- 
c+id). We obtain, taking into account the relations between the residues and the results of the sum 
rules shown in the previous section: 



-1-00 



co^-co 



a + ib 



- + ■ 



-a + ib 



c + id -CO -c + id-Q) 



J 



. ,y.( X 2{fO-iY^Q 2iS, 



(22) 



Z)(-fy) D(-fy) 



and 
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-l-oo 



(0^-0) 



{a + ibXc + idy {a + ibX-c + id) 



,2 A 



V 



c + id-Q) 



l{(0-iY)S^ 2iS, .2 (3)/ \ 



D{-co) D{-(o) 



c + id-Q) 



D{-(o) D{-(o) 



(23) 



We can separate real and imaginary part of the preceding expressions and obtain two couples of 
K.K. relations in the usual form, the first one originating from (22): 



\m{x''\(0^=-—P\d(0, 



2 

la 
jt 



1 2 2 

Q co^—co 



Tt 



So Im 



col -co ^ 




Tt 



-(5,-j<So)Re 



— (5,-j<So)lm 

Tt 



^ 1 ^ 

D{-co) 

^ 1 ^ 

Di^-co) 



, (24) 



and the second one originating from (23): 



(0^Rq{z'^'\(0))=-P \dco, 

Tt 



collm(x^'\o)S) 2co 



2 2 

C0[ - CO 



Tt 



Iml 



2t / (3)/ x\ 2w „ r , (Wf Re(ar* '(6^1)) 2<y ^ „ 
iyMm(y<^'(ft>))= P \dco,^ — ^ — ^-^ + — ^2 Re 

Tt i cot - co Tt 



,i>(-«), 

/ 1 



Tt 



D{-co) 



^ 1 ^ 

^D{-co)^ 
f 1 



— (53-}<S2)lm 

Tt 



D{-co) 
(25) 



Explicitating S2 and S3 with the help of expressions (12) and (15) we obtain the following K.K. 
expression in terms of the lowest moments: 



2 2 

<y, - CO 



2t /' (3)/' ^^ nf J «1 ReGr (6^1)) 269/' 2o 

co" \yc\x^^^ {co))= P\dco^ — 2 2 + — K-^o - )Re 



f 1 



2 , 
+ — fyg^i Re 

Tt 



D{-o)) 



' 1 ^ 

/ 1 



Tt 



co^S^lm 
(26) 



It is interesting to note the perfect simmetry between the two K.K. relations (24) and (25) with 
respect to the terms originating from non-holomorphic contributions: it can be observed a clear 
correspondence between So and Si in the (24) with respectively S2 and S3 in the (25). We can 
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observe that all the terms originating from non-holomorphic contributions are linear with respect to 
the Si withi=0,l,2,3. 

The nonlinear K.K. relations we have here obtained are more general than the ones given 
before, and in each of them we have also terms which are proportional to the intensity of the probe 
beam and are directly related to the two lowest moments of the susceptibility So and Si. Our 
procedure is a direct one, and avoids the use of the maximum entropy hypothesis [37,38] adopted by 
Peiponen [30,39,40,41]. 

Our K.K. analyses, as given by (24) and (25), are analytically complete and overcome the 
difficulties explicitely shown in ref. [30] (see for instance fig. 3.6 at p. 60 and fig. 5.5 at p. 92). In 
comparison with the results of the maximum entropy method shown in fig. 5.5 at p. 92 in ref. [30], 
our results reproduce exactly the relations between the real and the imaginary part of the 
susceptibility, since no hj^otheses or approximations have been made. 

V. CONCLUSIONS 

We wish to summarize the main results obtained above as follows. 
The anharmonic oscillator model has been used to derive general expressions for the nonlinear 
susceptibility with Pump and Probe optical beams. The lowest order susceptibility is 
shown to contain terms proportional to IE2P and also terms proportional to lEiP . The former are 
holomorphic in the upper complex plane with respect to the variable cOj and consequently obey sum 
rules and Kramers-Kronig relations of the type derived previously. The latter contain poles in the 
upper complex co, plane, and give contributions to the nonlinear sum rules and K.K. relations which 
had been formerly neglected. These contributions are shown to depend on the zero and first moment 
of the susceptibility only. 

The above contributions to the sum rules have been shown to depend on the derivatives of 
the potential at the equilibrium position, the third order susceptibility depending on the fourth order 
derivative and on the square of the third order derivative. 

The above results are expected to be of general significance, in spite of the simplicity of the 
model adopted, because they are only determined by the analytic properties of the response function 
due to time causality. It can be shown in fact that the leading holomorphic term here considered, 
Z^^^{co^',co^, -(0^,0)2), gives the same result in the most general physical system, treated with the 
Kubo response function [42], provided one substitutes the expectation values of the derivatives for 
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the values at the equilibrium position., and one considers the corresponding tensor which obtains 
from the spatial directions of the applied fields. 

The results here presented can be of interest for a detailed analysis of experimental data 
because they show explicitely the dependence on the intensity of the probe beam itself which must 
appear in the third order susceptibility z'^^^i^^i)- Experimentally it is also possible to separate the 
contribution proportional to the pump intensity by using modulation techniques so that the 
holomorphic and the nonholomorphic contributions can be separately measured. We suggest 
experiments where both the probe frequency and the probe intensity can be varied, so as to verify 
the appearence of the additional contributions in the sum rules and K.K. relations. 



ACKNOLEDGMENTS 

We are grateful to S. Scandolo and to G. La Rocca for suggestions and useful discussions. 



APPENDIX A 

In this appendix we present the general formula for the n'" order Polarization at the generic 
frequency rtOi+scOi, as derived from the iterative procedure described in ref.[33]: 



1 



t,p,q 



1 



X 



V(x) 



dx 



I tj \ 



nn^ — - 

j=l h=l J^h ■ 



(Ala) 



with the conditions: 
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l<j<l ,\<h,h'<tj 

(l < ij^h < ij+\,h' < and {h^h'^ [pj^/, ^ pjj,. ) or [qj ^, ^ qjj^, )) 



(Alb) 



The frequency considered (rcOj+sco^) imposes the following constraint: 



J=l h=\ 



(Ale) 



and the order to which Polarization is considered gives: 



/ 'j 
j=\ h=i 



(Aid) 



In order to consider Pump and Probe experiments we specialize the expression for the probe 
frequency imposing r=l and s=0 in expression (4) and define the following quantity, which 
immediately results to be limited to odd n : 



(A2) 



The above expressions contain terms which are predominant for an intense pump beam and a 
weak probe, and are holomorphic in the upper complex cOj plane. This corresponds to the 
contribution considered in the Kubo [42] expression and used in [26] to derive K.K. relations. Such 
terms are: 



iuiA EE*;.* 



-X 
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X 



V(x) 



otc 



{ I 'J 

yj=\h=\ 



nn 



i=\ h=l ' 



(A3a) 



with the conditions: 



\<j<l ,\<h,h'<tj 

(l < ij^h < ij+i,h' < and {h^h'^ [qj^f, It 



(A3b) 



and the usual constraints regarding rispectively order of Polarization and frequency 



/ 0- 
j=l h=l 



(A3c) 



J=l h=l 



(A3d) 



with an additional condition to select those terms in the summation (6a) which are not functions of 



CO, 



= 



(A3e) 



The main properties of above described terms are the following: 

i) they are proportional to the intensity of the pump beam to the highest possible power, the 
(«-l)/2*; 

ii) they are olomorphic with respect to the variable Wj in the upper complex half-plane; 

iii) they have the slowest possible decrease at infinity with respect to cOj (oc (O^ ) 
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The first and third properties give them the predominance over the other terms of the same order; 
the second and the third one allow us to write two K.K. relations because of Titmarsch's theorem, 

one for Z^"\o)i) and the second one for (OiZ^"^{<^i)^ shown in ref. [26]. It is also possible to 

prove that the summation in (6a) for every n results to be real by using the relation 

, which holds for every m, and rearranging the terms in thr summation so 

that conjugate terms appear in couple. 

The preceding considerations and the analysis of the asymptotic behaviour as obtained from 
K.K. relations and from expression (A3) allow us to derive sum rules for the moments of the 

Z^"''{cOi) up to the third moment of the imaginary part, which is the only one different from zero. 

Here we give its explicit expression: 




with the same constraints (A3b-A3c-A3d-A3e) in the summations. It can easily verified that the 
above results are in agreement with the general results obtained previously for the third order 
susceptibility [26]. 
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